Integer filling metal insulator transitions in the degenerate Hubbard model 
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We obtain exact numerical solutions of the degenerate Hubbard model in the limit of large di- 
mensions (or large lattice connectivity). Successive Mott-Hubbard metal insulator transitions at 
integer fillings occur at intermediate values of the interaction and low enough temperature in the 
paramagnetic phase. The results are relevant for transition metal oxides with partially filled narrow 
degenerate bands. 



71.27.+a, 71.30. +h, 71.20.Ad 

The understanding of strongly correlated electron sys- 
tems is one of the current challenges in condensed mat- 
ter physics. In recent years, this problem has received a 
great deal of attention from theorists and experimental- 
ist alike. Yet, our knowledge of even some of the basic 
features of the proposed model Hamiltonians remains to 
a large extent only partial, except, perhaps, in the one- 
dimensional case. In consequence, the interpretation of 
the experimental data of strongly correlated electron sys- 
tems has to remain only speculative in most of the cases. 
In regard of this situation, exact results on properties of 
model Hamiltonians in well defined limits is very desir- 
able. The relevant role played by the band degeneracy 
in models of strongly correlated electrons has been long 
and largely recognized Jl|-Q|. However, the systematic 
treatment of such models poses in general even greater 
technical difficulties than non-degenerate ones. This par- 
ticularly applies to numerical approaches which have to 
deal with the exponential grow of the Hilbert space. 

The goal of this paper is to demonstrate the existence 
of successive metal insulator transitions (MIT) at the in- 
teger fillings n — 1, 2, 3 in the two band degenerate Hub- 
bard model within the "Local impurity self-consistent 
approximation" (LISA) || which is exact in the limit of 
large dimensions (or large lattice connectivity) || . These 
MIT occur within the paramagnetic phase for intermedi- 
ate values of the interaction. We obtain exact numerical 
solutions of the model using an extension of the Hirsch 
and Fye quantum Monte Carlo (HFQMC) algorithm [@]|] 
for the solution of the associated impurity problem within 
the LISA. This impurity problem is a generalized single 
impurity Anderson model where the impurity and con- 
duction band operators carry an orbital index. 

The two band degenerate Hubbard model reads, 
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(ij) labels nearest neighbor sites and a, b = 1, 2 is the 
orbital index. This Hamiltonian is rotational invariant 
in spin and real space and the usual approximation U 
and J independent of band indices is assumed The 
parameter U is due to on-site Coulomb repulsion and 
the exchange parameter J originates the Hund's cou- 
pling. For simplicity, we shall further assume tfj = — tS a b 
and neglect the last "spin flip" term in (|l|). The re- 
sulting model Hamiltonian is relevant for for electronic 
systems with partially filled narrow degenerate bands. 
Examples of such systems are the 3d transition metal 
oxides Ri- X A X M O3 with three-dimensional perovskite- 
type structure, where R = La, Y, A = Ca, Sr and 
the transition metal M = Ti, V, Cr. Fujimori has con- 
structed a phase diagram by classifying these compounds 
as correlated metals or Mott insulators according to their 
electronic properties [ |To[ . Another and particularly no- 
table example is the extensively investigated V2O3 com- 
pound which has a MIT as a function of temperature, 
pressure and chemical substitution It is important 
to note that in bipartite lattices such as the hypercubic 
or the Bethe lattice the groundstate of the model has 
long range order B. In the single band case it is a spin 
antiferromagnetic (Neel) state, while in the degenerate 
model (with J w 0) one can have either a spin ordered 
and/or an orbital ordered state. Other types of orderings 
are also a priori possible, as for instance, ferromagnetic 
(for large J). The stability of the different phases will de- 
pend on the values of the various interaction parameters 
appearing in the Hamiltonian (Q). Since the parameter 
space is rather large, we defer the investigation of the full 
phase diagram for a later publication. Here we choose to 
concentrate on the interesting question of the destruction 
of a metallic state by the sole effect of electronic correla- 
tions (Mott transition). Therefore, we shall restrict our- 
selves to the paramagnetic phase (in both spin and orbit 
indices). It is important to note, however, that the para- 
magnetic solutions that we obtain are indeed the true 
groundstate in models that include next nearest neighbor 
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hopping (that lifts the rather artificial nesting property) 
and also in models with hopping disorder |12|T3l. 



In this paper we shall consider the model on a Bethe 
lattice of connectivity z — > oo which renders a semicircu- 
lar bare density of states p°{e) — (2/nD)y/l — (e/D) 2 , 
with D = 2t. Our choice is motivated from the fact 
that p°(e) shares common features with a realistic three 
dimensional cubic tight binding model, namely, it is 
bounded and has the same square root behavior at the 
edges. We shall set the half-bandwidth D = 1 and J=0 
(the case J ^ will be considered elsewhere). 

Within LISA, the lattice model is exactly mapped onto 
its associated impurity problem supplemented with a self- 
consistency condition This approximation becomes 
exact in the limit of the dimensionality d — > oo (or con- 
nectivity z — * oo). The resulting effective local action at 
a particular (any) lattice site reads, fl4|,|5|| 
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and the Weiss field is related to the local Green's 
function by the self-consistency condition, 



Q a aa (ioj) = iio + fi - t 2 G aa {iuj). 
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The LISA equations (|j)-(||) have to be self-consistently 
solved for. In practice this is done by iteration. One be- 
gins with a guess for Q® a which is used in (||) to solve the 
many-body problem and calculate G aa ■ This is then used 
as input in (||) to get a new Q^ a and the process is iter- 
ated until self-consistency is attained. Eq.|| also defines 
an impurity problem where the local site is coupled to an 
effective "conduction band" with hybridization function 
A(iu>) = t 2 G aa -(iuj) similarly to the Anderson single im- 
purity problem. This is the associated impurity problem 
within the LISA which has to be solved anew at each 
step of the iteration. To solve this degenerate impurity 
problem we use an extension of the HFQMC algorithm. 
The details for the implementation can be found in Ref. 
H . One of its salient features is the absence of the neg- 
ative sign problem which allows the investigation of low 
temperatures (~ 10~ 2 of the bandwidth). Converged 
solutions are typically obtained after 8 iterations using 
L = 64 time slices in a few hours on a workstation. 

We now turn to the discussion of the results from the 
solution of the LISA equations (§)-(§). In Fig.|l| we show 
the total occupation number n = (n^ + + ni\ + 
ni\) an d ano - it s first derivative k = dn/dfi (which is 
proportional to the compressibility), as a function of the 
chemical potential fj,. 
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FIG. 1. Occupation number n and k = dn/dfi (dotted and 
solid curves) as a function of the chemical potential fi. k is 
obtained by numerical differentiation. The data are obtained 
at U = 2.5, 3, 4 and T = 1/8 (top to bottom). Note: for better 
comparison u is measured respect the particle-hole symmetry 
point at | U. 

For U = 2.5 we observe drops in n at values of fj, cor- 
responding to the integer fillings n = 1,2,3 that indi- 
cates the onset of a correlated metallic state. This is 
demonstrated in Fig.|| where we plot the quasiparticle 
residue Z = (1 — 9ImE / duin)" 1 and the specific heat 
7/70 (= m* /mo in the d = 00 limit). We observe that as 
the system approaches the integer fillings the quasiparti- 
cles at the Fermi level become heavier as a consequence 
of the proximity to the MIT. On the other hand, at n = 
and 4 the quasiparticle residue approaches unity as the 
correlations vanish. It is worth noting the asymmetry in 
the behavior around n — 1 and 3. As one approaches 
n = 1 from below the specific heat diverges as the in- 
verse of the doping <5 _1 {8=1 — n) [^5). However, as 
we approach n = 1 from above the enhancement of 7/70 
seems to be faster. The results for < n < 1 are sim- 
ilar as the obtained within a one band model and gives 
a posteriori justification for its correct prediction of the 
behavior of 7/70 in La\- x Sr x TiOj, |Q. It would be in- 
teresting to test the asymmetry prediction of our model 
by approaching n = 1 from above in oxigen deficient 
LaTiOz-y- Around n = 2, which is the particle-hole 
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symmetry point, the enhancement of 7/70 is symmetric. 




FIG. 2. Quasiparticle residue Z and specific heat 7/70 
(bottom and top) as a function of the occupation number n 
obtained for U = 2.5 and T = 1/8. We estimate the slope in 
the self-energy E as the ratio of E/w n for the first Matsubara 
frequency. The wiggles are due to the QMC noise. 

Back to Fig.|l], we increase the value of the interac- 
tion strength to U = 3, and observe that the compress- 
ibility at 71 = 1 and n — 3 vanishes as the system 
goes through a metal insulator transition, [p~5f in fact, 
U c {n = 1) = U c (n = 3) due to particle-hole symmetry. 
On the other hand, the compressibility at n = 2 also de- 
creases but still remains finite, and the MIT is reached 
upon further increasing the interaction. Note that k also 
goes to zero at the endpoints \i — | « ±2U, but these 
are band insulating states which correspond to the com- 
pletely empty n = and completely full n = 4 state. 

We find that the MIT of the degenerate model at 
n = 1,2,3 have many features in common with the tran- 
sition that occurs in the single band case j^] . In partic- 
ular, at low enough temperature T < 0.04 for n = 1,3 
and T < 0.06 for n — 2, and within a finite interval of 
the interaction U two distinct solutions coexist: one has 
metallic character while the other shows a gap at low fre- 
quencies. In Fig.|] we show for the two fillings n = 1 and 
2 the coexistent Green's functions obtained at T = 1/32. 
To select either solution within the coexistent region, one 
has to choose an appropriate "seed" at the beginning of 
the iteration. These Green's functions are obtained as a 
function of Matsubara frequencies, i.e., they live on the 
imaginary axis. In order to obtain the density of states 
they have to be analytically continued to the real axis. 
Nevertheless, the metallic or insulating nature of the so- 
lution can be unambiguously determined from their low 
frequency behavior. In fact, the metallic Green's function 
behaves as ImG(w) 7^ for to — > and can be continu- 
ously connected to the non-interacting solution. On the 
other hand, the insulating Green's function shows linear 
behavior at small frequencies which is a signature of a 



gap opening in the density of states; thus, in this case, 
ImG(cij) — > for uj — > and can be continuously con- 
nected to the atomic limit solution (t — ► 0). 
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FIG. 3. Top: Imaginary part of the Green's function as a 
function of Matsubara frequency at T = 1/32, U = 2.9 and 
n = 1. The solid dots correspond to the metallic solution 
and the open dots to the insulating one. Bottom: Idem for 
T = 1/32, U = 3.4 and n = 2. 

Calling U c i(n) the minimum value of the interaction 
for which an insulating solution is allowed at given n, and 
U C 2(n) the maximum value for which a metallic solution 
is allowed; [|l2| we have established that at the lowest 
temperature investigated T — 1/32, Z7 c i(l) = U c i(3) w 
2.8, LT c2 (l) = t7 c2 (3) « 3 and U cl (2) « 3.1, U c2 (2) » 3.7 
Using similar arguments as for the single band model, |l6| ] 
one can demonstrate that at T = the metallic solution 
is lower in energy, thus, the T — MIT occurs at U C 2 and 
is of second order. A direct consequence of two solutions 
being allowed in regions of the (U, T) plane at given inte- 
ger fillings, is the existence of a first order transition line 
defined by the crossing of their free energies. This line 
starts at U C 2 at T = and has a negative slope due to 
the large entropy of the insulating state, while at finite 
T it ends at another second order point where the two 
solutions eventually merge. These results are condensed 
into the phase diagram of Figj| which shows qualitative 
agreement with that of Ref . |f(| . 
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FIG. 4. Interaction U vs. filling n phase diagram for 
the paramagnetic solution. The I s * order transition lines 
at n = 1,2,3 (dotted lines) end at 2 nd order critical 
points. Diamonds indicate the approximate position of the 
T — critical points (U C 2), and squares approximately in- 
dicate the finite T critical point where the two solutions 
merge (T(n = 1,3) « 0.04, T(n = 2) « 0.06). Solid 
lines at n = 1,2,3 indicate the Mott insulator states, and 
dashed-dotted lines at n — and 4 indicate the empty and 
full band insulator states. 

From the experimental standpoint, it has been argued 
by Fujimori et al. |lj] that the systematic analysis of 
the electronic properties of various transition metal oxide 
compounds with nominally one electron in a degenerate 
d— band shows evidence for a MIT due to electron cor- 
relations. Moreover, the extensively investigated V2O3 
compound contains nominally two electrons in a dou- 
bly degenerate e g band and its phase diagram displays 
a paramagnetic MIT line which is first order jflj . Never- 
theless, Castellani et al. |18 have argued that a relevant 
model for this system should contain just one electron in 
a doubly degenerate e g band. In either case the qualita- 
tive physics is being correctly predicted by our results. 

In conclusion, we obtain exact numerical solutions of 
the LISA mean field equations for the degenerate Hub- 
bard model. These are exact solutions of the model in 
the limit of d = 00. We find that within the param- 
agnetic phase the model has metal insulator transitions 
at the integer fillings n = 1,2,3 which are similar in 
character as the single band case. In particular we find 
that U c (n = 2) > U c (n = 1) = U c (n = 3), and that 
there are regions in the (£7, T) plane where two different 
solutions are allowed, one metallic and one insulating, 
which leads to a first order MIT line at finite tempera- 
tures. An important remark is that the observation of a 
similar scenario for the integer filling MIT in both, the 
single and degenerate band model, validates the often 
made assumption about the relevance of a single band 
Hubbard model for the qualitative investigation of com- 
pounds with narrow degenerate bands. This remark can 
be considered an a posteriori justification for the recent 



success in the interpretation of low frequency spectro- 
scopies of correlated systems with band degeneracy using 
a single band model [ pj|pO| ]. However, the previous state- 
ment has to be taken with care as we have seen that some 
quantities, as for instance the specific heat, may show a 
different behavior depending how the MIT is approached. 

Many interesting questions are open ahead, for in- 
stance, whether the present scenario is modified by the 
introduction of a finite J and the study of phases with 
long range order as function of U, J, T and n. These 
are certainly relevant questions for the understanding of 
the physics of transition metal oxides with partially filled 
d— bands. 

Valuable conversations with A. Georges, I. H. Inoue, 
H. Kajueter and G. Kotliar are gratefully acknowledged. 
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